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ABSTRACT 

The study of Lya emission in the high-redshift universe is a useful probe of the 
epoch of reionization, as the Lya line should be attenuated by the intergalactic 
medium (IGM) at low to moderate neutral hydrogen fractions. Here we present 
the results of a deep and wide imaging search for Lya emitters in the COSMOS 
field. We have used two ultra-narrowband filters (filter width of ~8 -9 A) on the 
NEWFIRM camera, installed on the Mayall 4m telescope at Kitt Peak National 
Observatory, in order to isolate Lya emitters at z = 7.7; such ultra-narrowband 
imaging searches have proved to be excellent at detecting Lya emitters. We found 
5cr detections of four candidate Lya emitters in a survey volume of 2.8 x 10 4 Mpc 3 
(total survey area ~760 arcmin 2 ). Each candidate has a line flux greater than 
8 x 1CU 18 erg s™ 1 cm" 2 . Using these results to construct a luminosity function 
and comparing to previously established Lya luminosity functions at z = 5.7 and 
z = 6.5, we find no conclusive evidence for evolution of the luminosity function 
between z = 5.7 and z = 7.7. Statistical Monte Carlo simulations suggest that 
half of these candidates are real z = 7.7 targets, and spectroscopic followup will 
be required to verify the redshift of these candidates. However, our results are 
consistent with no strong evolution in the neutral hydrogen fraction of the IGM 
between z = 5.7 and z = 7.7, even if only one or two of the z = 7.7 candidates 
are spectroscopically confirmed. 

Subject headings: dark ages, reionization — galaxies: high-redshift — galaxies: 
luminosity function 
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INTRODUCTION 



Direct observations of distant galaxies remain the most straightforward way to probe 
the fundamental nature of the high-redshift universe. Such observations can provide some 
much-needed constraints on numerical simulations, which may provide better answers to the 
question of how large scale structure forms and how star formation begins in dark matter 
halos. Star formation in early galaxies is dependent upon the mechanism by which gas 
cools; that mechanism itself is dependent upon the ionization state and metal enrichment 
of that gas which are not well constrained at high redshift. The Lya emission line is a very 
useful tool for the detection of high-redshift galaxies, as the earliest stars in the universe 
sho uld ionize surround ing hydrogen gas, which will then recombine to produce Lya emission 
(see lWillis et al.l 120081 ). This Lya line can be probed quite effectively at high redshifts via 



the use of narrowband filters, w hich focus on regi ons with low sky background and that 



are free of strong OH lines (e.g., ICuby et al.l 120071 ). High-redshift objects should have es- 



sentially no flux blueward of rest-frame Lya and none blueward of rest-frame 912 A; this 
is a result of th e Lya forest effect at high redshift, due to strong absorption by interven 
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20091 ) and have resulted in samples of galaxies over a range of redshifts, includ ing the spec- 
troscopic confirmation of a Lya emitting galaxy at z — 6.96 (live et al.l 120061 ). Even when 
these narrowband surveys do not successfully detect high- z objects, su c h null results can 



be used to constra in the Lya luminosity function (e.g., ICuby et al.l 120071 ; IWillis et al.l 12008 
Sobral et al-lbopgh . 



The early universe is expected to be metal poor, but meta ls have been detected in the 
intergalactic medium (IGM) at z — 5.7 (IRyan- Weber et al.ll2006l ). and so the IGM must have 
been enriched in metals by z ~ 6 at the latest, wit h recent results tent atively indicating a 
metallicity downturn between z ~ 5.7 and z ~ 5 (jSimcoe et al.l 120 111 ). Additionally, the 
Lya line is sensitive to IGM obscu ration at neut r al hydrogen fractions ranging from low 
to high (e.g., 10% < x HI < 100%; iHaimanl |2002| ; ISantosI 120041 ). and thus observations of 
Lya -emitting galaxies serve as a powerful probe of the reionization history of the universe; 
the Gunn-Peterso n test, for example, is only useful when the neutral gas fraction is <1% 
( ICuby et al.l 120071 ). Increasing the neutral hyd rogen fraction in t he IGM increases the at- 
tenuation of Lya emission from those galaxies (jSobral et al.l 120091 ). As this neutral fraction 
increases, the Lya luminosity function will vary according to the amount of light being at- 
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tenuated by the IGM. Previous estimates of the redshift of the epoch of reionization using 
constraints from studies of Lya emitters do not concur with constraints derived from polar- 
ization observations of the CMB. The latter suggests that the redshift of reionizatio n is z re = 
10.5 ± 1.2, should reionization be an instantaneous process ( IKomatsu et al.ll201l[ ). whereas 
the former have indicated a significantly later end to the epoch. Constraints on the Lya lu- 
minosity function (LF) can assist in determination of the redshift at which reionization has 
been completed, owing to the resonant scattering of Lya photons in a neutral IGM. If the 
intrinsic number density of young galaxies remains constant over redshift, then a significant 
decline in the observed Lya LF at a given redshift could indicate a change in IGM phase. On 
the lower redshift end, Malhotra & Rhoads (2004) found no significant evolution of Lya LF 
between 5.7 < z < 6.6; whereas at higher reds hifts, an evolution of Lya LF b etween 6.5 
< z < 7 is suggested based on single detections (live et al.ll2006l ; lOta et al.ll2008l ). Ouchi et 
al. (2008) found little evolution between z ~ 3 and z ~ 6 in the observed LFs, although 
they suggest a real evolution, with increase in intrinsic Lya luminosity being canceled out by 
increase in IGM absorption. Curtis-Lake et al. (2011) have recently identified Lya emitters 
at a high rate over 6 < z < 6.5 in UKIDSS. Ono et al. (2011) have measured a decrease in 
Lya emission line detection fraction over 6 < z < 7, as have Schenker et al. (2011) over 6 < 
z < 8 and Pentericci et al. (2011; following up work by Fontana et al. (2010)) over 6 < z < 
7; these studies are UV-continuum-selected galaxies (i.e., Lyman break galaxies) rather than 
selected via Lya emission line but share the same goal, and all conclude that the neutral 
hydrogen fraction of the IGM is increasing over those epochs. The seven Lya candidates 
found by Hibon et al. (2010) at z — 7.7 would indicate no evolution of the Lya LF if they are 
found to be spectroscopically confirmed (although soon to be published work suggests that 
at least five are not confirmed; see Clement et al. (2011)); the four Lya candidates found by 
Tilvi et al. (2010) could either indicate mild or no evolution between 6.5 < z < 7.7 depend- 
ing on the number of candidates which are confirmed. Two Lyman break galaxies at z >7.0 
were recently spectroscopically confirmed by Vanzella et al. (2011), their luminosities fairly 
consistent with that of the well-known confirmed z = 6.96 lye et al. (2006) Lya emitter. At 
present, small number statistics severely affect our ability to draw definite conclusions on 
the Lya luminosity function and properties of the IGM. It is therefore essential to expand 
the sample of high-redshift candidates at epochs when the IGM became metal-enriched and 
reionized, and thus shed light on the nature of the early universe. 

As redshift increases, galaxy sizes and luminosities decrease, and cosmological dimming 
must als o be taken into account. Because of this, detection of galaxies at z > 7 can be quite 
difficult ([Ferguson et al.l 12004 ; iBouwens et all 120061 ; ICapak et al.ll201ll ). Thus it is essential 
to have a large survey volume in order to detect a sufficient number of high- z objects. This 
conclusion is bolstered by the biased nature of galaxy formation and the non-uniformity of 
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have significant depth but small area (e.g., iBouwens et al.ll2010t lOesch et al.l 2010) o r have 



small volume but high magnification by virtue of cluster lensing ([Richard et al 



2007, 2008 



Stark et al.ll2007l ). so all results are affected by cosmic variance. There is a need for surveys 
which probe both a deep and wide region in order to best constrain global properties of 
high-z galaxies; the present paper reports on a survey that tries to fill that role. 

In this paper, we present the results from a search for Lya emitting galaxies at z = 
7.7 in the COSMOS field, utilizing custom-made ultra-narrowband filters which are tuned 
to avoid the OH sky lines and thus reach extremely low infrared sky backgrounds. The 
organization of this paper is as follows. In Section [21 we discuss our observations and 
reduction of the data that we obtained, as well as photometric calibration. In Section [21 
we describe our method of candidate Lya emitter selection and the basic properties of the 
resulting candidate Lya emitters. In Section HJ we discuss possible sources of contamination 
in our samples. In Section [5j we estimate the number of Lya emitting galaxies that we 
should expect to find in our survey through the use of a detailed Monte Carlo simulation. In 
Section [HI we present the Lya luminosity function derived from our candidates and compare 
to previously derived Lya LFs. Finally, in Section [TJ we summarize our conclusions. This 
work shares authors, instrument, and technique with the work of Tilvi et al. (2010), and 
can thus be viewed as part of a series with that paper. Throughout this work, we assume 
a standard flat ACDM cosmology with Q m = 0.3, Q\ = 0.7, and h = 0.71, where Q m , Q\, 
and h are the matter density, dark energy density, and Hubble parameter (in units of 100 
km s -1 Mpc -1 ), respectively. All magnitudes listed are in the AB magnitude system. 



2. OBSERVATIONS AND DATA 

2.1. Observations with NEWFIRM 

Our observations were centered on the Cosmological Evolution Survey (COSMOS) field 
(RA 10:00:28.6, Dec. +02:12:21.0), taking advantage of the large amount of ancillary data 
available on this fielcfl. We used the NOAO Extremely Wide-Field Infrared Mosaic (NEW- 
FIRM) camera on the Mayall 4m telescope at Kitt Peak National Observatory (KPNO) on 
three different observing runs over the course of two years (2008 February 28 - March 14; 
2009 January 29 - February 1; 2009 February 17 - March 2) for a total of roughly 100 hours 
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over 32 nights. Average seeing over the course of these exposures was ~1.2". 

The NEWFIRM camera, sensitive to 1-2.4 /iin wavelengths, is a wide- field imager con- 
sisting of four mosaiced 2048 x 2048 pixel ALADDIN InSb arrays, 0.4" pixel -1 , for a cu- 
mulative field of view of 27.6' x 27.6' (cumulative area ~760 arcmin 2 ). Our observing time 
was split evenly between two University of Maryland custom-made ultra-narrowband (UNB) 
filters (R~1000) centered at 1.056 and 1.063 /mi, with FWHM of 7.4 and 8.1 A, respectively. 
These UNB filters were designed in order to isolate Lya emitters at z ~ 7.7, while simul- 
taneously avoiding the bright OH lines near these wavelengths (Figure [T] shows UNB filter 
profiles). The transmitted wavelength varies across the field according to roughly mX = 
fid cos 9 (here the order, m, index of refraction, fi, and thickness of the filter, d, are constant; 
9 is the angle of incidence of light onto the detector); the path length increases as one gets 
further off-axis and the wavelength of transmitted light is correspondingly shorter. This 
effect is more apparent for a UNB filter in a wide field imager, since the bandpass is narrow 
and the range of angle of incidence is large (see Figured]). Atmospheric absorption (primarily 
due to O2 and H 2 in the IR) is irrelevant in these bandpasses. Each 1200s exposure science 
frame was taken using Fowler 8 sampling and we utilized random dithering within a 45" box 
after each exposure. 

Additional ancillary data in broadband filters were required for Lya candidate selection. 
We made use of publicly-available data from the COSMOS archived for this purpose: B, r, 
and % band data from Subaru and J band data from UKIRT. 



2.2. Data Reduction 

The reduction of our data was done by the NEWFIRM Science Pipeline (see Dickinson 
and Valdes 2009; Swaters et al. 2011). Frames with seeing of > 1.5" were rejected outright. 
The pipeline flagged pixels affected by detector blemishes, saturation, and persistence, and 
then subtracted the dark current, linearized the data, and applied the dome flat. Image 
gradients remaining after the flat fielding were subtracted out. The astrometric solution was 
determined from 2MASS stars in the field; these same stars were also used to determine 
an initial photometric calibration of the data. In the pipeline, the sky was subtracted in 
a two-pass approach. First, the sky was subtracted using a running- median window, and 
then the data were combined by taking the median over all the exposures. This so-called 
harsh stack was then used to identify and mask sources. The mask was then applied to the 
original images, and the sky-subtraction was repeated. Cosmic ray hits and other transient 
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phenomena were detected by comparing individual images against the first-pass stack, and 
outliers were flagged. 

Starting with these pipeline products, custom IDL4 scripts designed by Krug & Swaters 
were used to eliminate artifacts such as OH rings and striping due to data readout from 
the science frames. For OH ring elimination, pixel values were separated into 1000 radial 
bins across one NEWFIRM chip at a time. Owing to the unevenness of the NEWFIRM 
chip gaps, the pixel center of the OH ring varied from chip to chip; the ring centers were 
determined for each chip through a combination of visual inspection and fitting rings to the 
values at each pixel. Data within the radial bins were then median smoothed, and these 
smoothed bins were then subtracted from the original image data. Stripe removal was also 
performed on individual NEWFIRM chips at one time. Each chip was divided into horizontal 
or vertical strips, depending on data readout direction. Pixels across each strip were summed 
and averaged, and the original data pixels in each strip were subtracted by the average pixel 
value in that strip. A comparison of stacks made from frames before and after ring and 
stripe removal showed that the signal-to-noise ratio (SNR) increased by 1.2 times on average 
across the frames, with SNR increasing by as much as 4-7 times in regions where the OH 
rings were strong prior to algorithm implementation. The SNR in the most extreme rings, 
however (namely those areas as seen in Figure [1] where the OH lines cross the outer edge 
filters, particularly the strongest OH ring in the 1.063 /im filter), was not quite high enough 
for candidate detection (see Section [3]). 

Following removal of artifacts and sky subtraction, the world coordinate system of each 
individual frame was assigned using the IRAF task n fives, which matches as many sources 
as possible to those of the 2MASS point source catalog, searching within a 450" radius. The 
mscimage task was used to reproject all four chips onto a single image with the same pixel 
grid. Finally, images were combined into two yearly stacks (weighted by seeing) as well as 
one full data stack for each filter via imcombine; these were median combined using noise 
values from mscstat and zeropoint magnitudes as calculated based on the 2MASS catalog 
during the prior nfwes task. This resulted in two full UNB stacks - one each for the 1.056 
and 1.063 /im filters - as well as stacks from year one and year two for each filter (in order 
to properly check for transients). Broadband data mosaics supplied by the COSMOS team 
were already reduced. 



2.3. Photometric Calibration 



Photometric calib ration was performed by m aking use of the 2MASS point source cat- 
alog. SExtractor (SE; iBertin fc Arnouts I (|l996h ) was run on each NEWFIRM UNB stack 
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using a magnitude zero point set to zero. The resulting catalogs were matched with the 
2MASS point source catalog; in order to avoid saturated sources or objects that are too 
faint, all stars with 2MASS J band magnitudes (AB) of less than 13.0 or greater than 16.5 
were discarded. Plots were made comparing the SE UNB magnitudes with 2MASS J mag- 
nitudes in order to determine an appropriate zeropoint. Before SE could be run with the 
proper magnitude zeropoint, it was necessary to determine a color correction between the 
2MASS J magnitude and the NEWFIRM UNB magnitudes, as they do not share the same 
central wavelength. 2MASS J, H, and K band magnitudes for overlapping stars were con- 
verted to flux and used to calculate the flux expected at 1.056 and 1.063 /mi, the central 
wavelengths of the NEWFIRM UNB filters. The difference between 2MASS J flux and ex- 
trapolated NEWFIRM UNB flux corresponds to a magnitude difference of ~+0.1 mag in J. 
This correction was then applied to the previously determined magnitude zeropoint. 



2.4. Limiting Magnitude 



In order to determine which of our detections can be considered real, we must obtain 
some estimate of the limiting magnitude of our science frames. We have chosen to define our 
limiting magnitude as the 50% completeness limit in our frames. To determine this, we inject 
artificial sources into our frames. These sources are randomly distributed throughout each 
frame, although care is taken to avoid locations within 4" of any existing bright star. Two 
hundred artificial point sources are generated in each of 0.1 magnitude bins, ranging between 
20 and 24 mag for UNB frames, 21 and 25 for J band, and 24 and 28 mag for the broadband 
optical frames. Once these sources were generated, SE was run and the resultant SE catalog 
was matched with the catalog of artificial sources to determine a recovery fraction. This 
procedure was then iterated 25 times for each band. Once all iterations were completed, 
histograms of recovery fraction of artificial sources were plotted for each band in order to 
determine the magnitude at which 50% of all artificial sources were detected by SE (Figure 
[5] for UNB recovery fraction histograms). We define the limiting magnitude for each band 
as the 50% completeness limit for that filter, as this is the magnitude brighter than which 
sources can be reliably detected in these fields via our methods. Limiting magnitudes for 
the two NEWFIRM UNB bands were determined to be 22.4 and 22.5 (AB) for 1.056 /an 
and 1.063 /mi, respectively. Shifts in central wavelength of the UNB filters with distance 
from the field center were not found to significantly affect these limiting magnitudes. For 
the broadband data, limiting magnitudes were found to be 27.3, 26.6, 26.3, and 24.0 for B, 
r, i, and J, respectively. 
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3. CANDIDATE Lya EMITTERS 

3.1. Selection Criteria 

To select potential candidate Lya emitters, we initially ran SE on the NEWFIRM UNB 
frames individually, with minimum detection area of 5 pixels and detection threshold of 
1.2 sigmas. We then used dual- image mode in SE, which takes the coordinates of objects 
detected in the NEWFIRM UNB frames but measures fluxes at those coordinates in the 
UKIRT J band image. Thus we obtained J band fluxes at the exact locations of NEWFIRM 
UNB targets. For all of our dual-image mode detections, we made use of the SExtractor 
auto aperture size scaling. We also needed to obtain optical fluxes at those coordinates, 
so in order to do so, we created a chi-square d image of a combi nation of B, r, and % band 



data from Subaru using the program Swarp (IBertin et al.ll2002l ). The chi-squared image is 



created by deriving the distribution of the pixels that are dominated by object flux rather 



than those dominated by sky noise (jSzalay et al.lll999l ). This is an optimal combination of 



images which enhances real objects in the resultant frame rather than a simple stack. We 
then ran SE in dual-image mode on this combined optical image in the same vein as the J 
images. We also ran dual-image mode on the two NEWFIRM UNB frames together - first 
using the coordinates for 1.056 /im detections but getting 1.063 fim flux, and vice versa - 
as any potential candidate Lya emitter should appear in one of the UNB frames but not 
both. We likely miss highly extended Lya blobs with this technique, especially following sky 
subtraction, but such highly extended sources would be difficult to isolate in any case. 

Once the SE catalogs were created, we ran a custom python script to comb through 
those catalogs and identify potential candidates. There were five main initial criteria that 
had to be met in order for an object to be considered a candidate (J data refers to the 
UKIRT band, centered at 1.2/im and outside the UNB filter range): 

1. UNB signal-to-noise ratio (SNR) of at least 5. 

2. UNB excess /„(UNB) - f v {J) of at least 3a. 

3. Flux ratio /„(UNB) / f u (J) of at least 2. 

4. UNB SNR in the other UNB filter of no more than 2 (e.g. if detected in 1.056 /mi, 
must not have a 1.063 fim SNR of more than 2). 

5. Optical SNR in the chi-squared broadband image of no more than 2. 



These selection criteria have been utilized in lower redshift (z = 4.5, 5.7) searches for 
Lya emitters and have detected those emitters at a 70-80% success rate after spectroscopic 
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confirmation, and we are confident that these criteria will translate to higher redshift, as 
the fundamental physics of the Lya forest should not < 
used in a successful spectroscopic search at z = 6.96 



e.g., 


Rhoads & Mall 


lotra 


2001; 


Dawson et al. 


2007 




Wane et al. 



2009h . These criteria are also the same as those used in the work of Tilvi et al. (2010). 



The first three criteria are used to isolate emission line sources. The remaining two are used 
to eliminate as many low-redshift sources as possible; as detailed in Section [fl Lya emitters 
should have no flux blueward of their Lya emission line, and the Lya line should be narrow 
enough that it is only detected in one NEWFIRM UNB filter and not both. 

Following execution of the selection script, we were left with 65 potential candidates out 
of 31254 initial detections for the 1.056 /im band, and 110 potential candidates out of 32382 
initial detections for the 1.063 /j,m band. We then matched our target lists to the existing 
COSMOS source catalogs (obtained from Peter Capak of the COSMOS team; personal 
communication). All targets which met our selection criteria but which were also in the 
COSMOS catalog with a measured photometric redshift were set aside (Section ETTl for further 
information on these low-redshift interlopers). Once the low- 2 interlopers were removed, 
target lists were then narrowed down further through a variety of methods: eliminating all 
targets on the edges of the chip which lie in or outside the main OH-line ring; eliminating 
all targets that lie within two arcseconds of the chip gap (5 pixels); eliminating all targets 
within two arcseconds from a very bright star. Roughly 90% of the initial non-interloper 
candidates were removed this way. Basic visual inspection was then performed as a sanity 
check. Our targets were also compared with the yearly stacks for each filter (both via SE 
and visual inspection) in order to ensure that these targets were not transients (Section HJ). 
As a final sanity check, we did not consider any candidates with magnitudes fainter than the 
limiting magnitudes for each filter (Section 12. 4p . 



3.2. Results 

After all these tests were completed, we were left with a total of 4 candidates brighter 
than the 50% completeness limit - 3 candidates in 1.056 /xm, 1 in 1.063 jttm. Three of these 
candidates lie at the survey line flux limit. The coordinates and basic properties of these 
four luminous candidates - AB magnitude, line flux, and luminosity - are listed in Table HJ 
Line flux was calculated from the SExtractor magnitude using: 

F = l -°- 4 *( ma SAB+48.60) C w ergs -l cm -2 s -l 5 (1) 
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where magAB is the magnitude from the isophotal fit as reported by SExtractor, c is the 
speed of light, A is the central wavelength of the given filter, and W is the filter width (6.95 
xlO -4 /im for the 1.056 /im filter, 7.49 xl0~ 4 /im for 1.063). Whereas we used the SE auto 
fit for initial dual-image mode detection, we use isophotal fit for flux calculation to avoid 
losing signal and gaining noise. We can also estimate the star formation rates for these 
objects, assuming that that there is little attenuation of the Lya line by the neutral IGM 
and that the dust content along the line of sight is lo w. We use the f ollowing prescription 
from Ota et al. (2010), which uses the Kennicutt law (jKennicuttlll998l ) and assumes case B 
recombination: 

SFR(Lya ) = 9.1 x 10" 43 L(Lya ) M yr -1 . (2) 

Using this calculation, we find that star formation rates for our four Lya candidates range 
from 5 to 7.6 M Q yr -1 . Postage-stamp images of these candidates are shown in Figure [3J 
Note that each candidate is visible in one particular NEWFIRM UNB stamp but is not 
visible in the other NEWFIRM UNB stamp. Additionally, no candidate is visible either in 
the chi-squared broadband optical stamp nor the UKIRT J band stamp. As stated in Section 
[31 in order to eliminate most foreground galaxies, it is required that none of the candidates 
are detected in the broadband image; non-detection in the J band image is simply due to the 
faint continuum of the Lya candidates. As can be seen in Figure [3j the PSF is considerably 
better in the Subaru band stack than in our UNB data. To try to determine whether this 
could affect our candidate selection, we have attempted a modified aperture correction using 
a comparison between total SE auto flux in the Subaru field and flux derived from isophotal 
fits; this correction allows us to account for faint sources in the optical image. Following 
this correction, we re-ran our selection criteria. Three of our four high-z Lya candidates 
still passed the selection criteria, but candidate #3 failed selection following this aperture 
correction (optical SNR > 2). We still include candidate #3 in the following analysis but 
flag it in Table [2] as more uncertain than the other candidates. 



3.3. Lyct Equivalent Widths 



As we have a sample of four strong candidate Lyct emitters, it is worthwhile to compare 
the equivalent widths (EWs) of their Lya emission lines to those already noted in the litera- 
ture. Several published studies have spectroscopically identified Lya emitters with rest-frame 
equivalent widths of EW rest > 240A at z — 4.5 a nd 5.7, significantly higher than predicted by 
theoretical simulations of star-forming galaxies 
20061 : bawson et al1l2007l : iGronwall et "aDl2007l : buchi et al.ll2008h . 
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The calculation for rest-frame EW for our Lya emitters makes use of the fluxes measured 
in both the NEWFIRM UNB filters and the broadband UKIRT J data, as the Lya line does 
not appear in the J band (again, the J band data is centered at 1.2/im, with no overlap 
with our UNB filters). However, none of our candidates were detected in the J band, and 
thus we must use the J band limiting magnitude (Section I2.4[) to calculate an upper limit 
on J band flux. In this case, given the limiting magnitude of 24.0 in the UKIRT J band, 
the upper limit on the J band continuum flux f\ j = 1.9 x 10~ 19 erg s _1 cm" " 2 A" 1 . This 
value is used to calculate a lower limit on the rest-frame Lya EW for our candidates: 



EW rcst = — — x — — , 3 
fx,j 1 + z 

where /unb is the UNB line flux in erg s -1 cm -2 . This calculation assumes that the Lya line 
falls entirely within the transmission profile of our UNB filters. In reality, the UNB filters 
may not enclose all of the Lya emission. Moreover, we are using upper limits on the J band 
fluxes. The EWs measured in this way should thus be viewed as l-a lower limits. This 
calculation also assumes an exact redshift of z = 7.7. 

The resultant lower limits on the rest-frame EW (Lya ) are 7.32, 5.17, 4.89, and 4.84 
A for candidates 1-4, respectively. Thus our lower limit on EW for our candidates is EW res t 
> 4.8 A, much more consistent with theoretical predictions than the numbers quoted above 
for lower redshift surveys. Our lower limit EW is smaller than that of Hibon et al. (2010) 
by several angstroms, primarily owing to the difference in bandwidths of UNB filters used 
by these two surveys. This line width is also larger by ~2 A than the results of the Tilvi 
et al. (2010) survey, although this difference can be accounted for by the increase in depth 
of our J band data (limiting magnitude of 24.0 for our data versus 23.5 for their data). As 
future surveys obtain deeper J band data, we should be able to better constrain the lower 
limits of Lya EWs for these emitters. 



4. POSSIBLE SAMPLE CONTAMINATION 

There are several possible sources of contamination in our sample of candidate Lya emit- 
ters. These include such real sources as foreground emitters, transients, and cool L & T 
dwarfs, as well as false sources such as detector noise spikes and general false detections. We 
discuss each class of contaminants below. 



4.1. Foreground Emission Line Sources 



There are three main species of foreground emission line objects which are most likely to 
contaminate our sample, as each species should have a strong emission line that falls within 
our UNB filter window and, assuming faint continuum emission, negligible flux in nearby 
blue- and redward bands. These species are Ha A6563 emitters at z — 0.62, [O III] A5007 
emitters at z = 1.12, and [O II] A3 72 7 emitters at z = 1.85. Given the extensive amount of 
ancillary data available for the COSMOS field, we were able to first check our catalog with 
the main COSMOS catalog and eliminate any \ow-z interlopers from our Lya candidate list, 
as mentioned in Sectional A total of 3 interlopers fulfilled the five main Lya emitter selection 
criteria (Section[3]), all in the 1.063 /mi filter - one brighter than the limiting magnitude of the 
field - which were removed from the sample of Lya candidates. These faint objects, however, 
are not at the exact redshifts listed above: the COSMOS photometric redshifts are z ~ 1.5 
± 0.2 (lcr; possibly corresponding to redshifted FTy) , z ~ 1.7 ± 0.2 (possibly [Oil] A3727), 
z ~ 2.5 ± 0.2 (possibly broad Mg II A2798), but the errors may be even larger. It is possible 
that there remain some foreground emission line objects in our Lya candidate sample that 
do not have tabulated photometric redshifts. Our survey benefits greatly from the multi- 
wavelength coverage and accurate photometric redshifts of the COSMOS data catalog; the 
presence in our sample of interlopers with photo-zs different from the three expected species 
listed above may imply that surveys in other fields without such comprehensive ancillary 
data coverage suffer from greater low- z contamination than estimated. 

To estimate the number of additional possible emission line source interlopers that could 
remain among our candidates, we must estimate the minimum equivalent width required for 
these emission lines to contaminate our sample. We must also use the depth of our UNB 
image in order to calculate the minimum luminosities of these emission lines. As described 
in Section [23J the limiting magnitude of our 1.056 /xm stack is 22.4, which is equivalent to 
a 50% completeness limit in flux of 7.4 x 10~ 18 ergs -1 cm~ 2 . Given the redshifts of the 
aforementioned emission lines and using CosmoCalqj to calculate luminosity distances, we 
find that the minimum luminosities required to detect these emitters are 1.15 x 10 40 erg 
s~\ 4.96 x 10 40 erg s -1 , and 1.72 x 10 41 erg s" 1 for Ha, [O III], and [O II] respectively. To 
calculate the necessary minimum equivalent width (observer- frame) , we use the prescription 
from Rhoads & Malhotra (2001): 

EW min = - l] AA nb = \p± - ll AA nb , (4) 



3 http : //www. astro .ucla.edu/ ^wright/CosmoCalc .html 
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where / n b and /bb refer to the flux densities of the UNB and the chi-squared optical image, 
respectively, <r n b and a"bb refer to the flux measurement uncertainties in the two frames, and 
AA n b is the width of the UNB filter. For these calculations, we have simply used the 1.056 
/iin UNB filter. It is safe to assume that the continuum contribution to the overall flux in 
the UNB filter is negligible, since the overall width (effective FWHM) of the transmission 
profile of the 1.056 jjm filter is only 9 A (6.95 A). For our filters, <r n b = 1.06 x 10~ 29 erg cm -2 
s _1 Hz -1 and Obb = 2.1 x 10 -31 erg cm -2 s _1 Hz -1 . This results in a minimum equivalent 
width for the foreground emission line contaminants of EW m i n > 870A. 

Although the distribution of equivalent widths of these three species of emitters has not 
been probed at the high redshifts of our observations, we can scale the published results under 
the assumption that the luminosity functions of these species have not evolved significantly 
between the relevant redshifts. One particular recent study has probed emission line sources 
of Ha at z = 27, [ O III] at z = 0.51, and [O II] at z = 1.0 in the GOODS-South field 



( iStraughn et al.ll2009l ) . The minimum EW calculated above was scaled appropriately in each 
instance to the EW one would expect at the redshifts probed by the Straughn et al. (2009) 
survey. Using that value and the required minimum luminosity listed above, we determined 
how many emission line objects we should expect in our field, after scaling appropriately by 
the ratio of survey volumes. We find that we expect less than one additional interloper for 
each species (0.1 each for Ha and [O II], 0.3 for [O III]). If we relax the minimum equivalent 
width criterion to 500 A prior to scaling EWs to what would be expected at z = 7.7, the 
expectation increases to 1 additional interloper each for Ha and [O II], and 2 additional 
interlopers for [O III]. 



4.2. Other Possible Contaminants 

The utilization of individual yearly stacks, as well as nightly stacks, which were produced 
by the data reduction pipeline, allowed us to eliminate both possible detections due to noise 
spikes as well as real sources such as transient objects. We required during our candidate 
selection process that potential candidate Lya emitters be detected in all stacks as well as 
the overall stack, though the detection requirements for individual nightly/yearly stacks were 
less stringent. This eliminates contamination by noise spikes from the detector, as such noise 
spikes should not be present at the same coordinates across multiple nights over a span of 
years. Transients will also be eliminated by this yearly stack check, as supernovae should 
only be visible for a few weeks, not years. We were able to eliminate upward of twenty 
contaminants via this requirement. 
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In order to determine whether any false detections were contaminating our candidate 
list, we checked to see whether taking the negative of our UNB stacks would result in a 
detection. We used IRAF to multiply each UNB stack by -1 and then ran SExtractor on 
the resultant negative images. The SE negative stack catalogs were then run through the 
exact same selection process as the positive stack catalogs, but no candidates were identified 
this way. We are thus confident that the probability of false detections contaminating our 
Lya sample is insignificant. 

The final source of potential contaminants consists of cool stars - namely, L and T dwarfs 
- passing our selection process . We use previousl y observed relationships between spectral 



type and absolute magnitude (ITinney et al.ll2003l ) to determine in what J band magnitude 



range L and T dwarfs fall, and then we use our calculated J band limiting magnitude to 
determine distance ranges at which we should be able to see these objects. We find that L 
dwarfs could be detected between roughly 550 and 1700 pc, and T dwarfs could be detected 
between roughly 200 and 750 pc. These L and T dwarfs are m ost likely found within 



a Galactic disk scale height of 350 pc, however (IRyan et all 120051 ) - the number density 
drops significantly above the scale height of the disk - and thus knowing that these dwarfs 
typically have a volume density of no more than ~3 x 10 -3 pc~ 3 , we thus can conclude that 
we should expect at most one L or T dwarf in our survey. If we then take our selection 
criteria into account, we can determine if any L and T dwarfs would satisfy our narrowband 
excess criterion. Tilvi et al. (2010) used existing observed spectra of L and T dwarfs and 
a calculation of expected flux through our NEWFIRM filters. They determined that the 
flux which would be transmitted through the NEWFIRM UNB filters, in comparison to that 
which would be transmitted through the NEWFIRM J filter, would not be sufficient to pass 
our narrowband excess criterion. Thus, we should expect that no L nor T dwarfs would pass 
our selection criteria and contaminate our sample. 



5. MONTE CARLO SIMULATIONS 

As discussed in Section HJ we do not expect all of our four candidate Lya emitters to be 
real, as it is possible that further low-z emission line source interlopers may pass our selection 
criteria. In order to determine how many of our Lya candidates we should expect to be real, 
we run Monte Carlo statistical simulations using previously known Lya luminosity functions 
at lower redshift. In this work, we make use of the Kashikawa et al. (2006) Lya luminosity 
function at z = 6.5 and assume that there has been no significant evolution between z = 6.5 
and z = 7.7. These simulations also make use of the fact that the NEWFIRM UNB filter 
transmission curves may not encompass the full width of the expected Lya emission lines 
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in these objects, and thus our measurements may in fact be underestimating the Lya line 
flux produced by our Lya candidates. The transmission curves of our two UNB filters were 
obtained from the NEWFIRM websita9@. 

To begin our Monte Carlo simulations, we utilized the aforementioned Kashikawa et al. 
(2006) luminosity function to generate one million random galaxies. These galaxies were 
distributed with both a random Lya luminosity in the range of 10 42 erg s _1 to 1.5 x 10 43 
erg s _1 , as well as a random redshift in the range probed by our filters. In the case of the 
1.056 /mi band, this range is 7.66 < z < 7.71; for the 1.063 /mi band, this range is 7.72 < z 
< 7.76; we chose these ranges to correspond to the full width of the transmission profile at 
which the transmission falls to 5%. 

From these luminosities and redshifts, we were then able to assign a flux to each of our 
randomly generated galaxies. We made use of an asymmetric Lya line flux profile based 
on spectra of z = 5.7 Lya emitters taken by Rhoads et al. (2003). We were able to then 
use the detailed NEWFIRM UNB filter transmission curve profiles to derive the flux that 
would pass through the filter via a convolution of filter profile and line flux according to 
/trans = / f\T\d\, where T x is the NEWFIRM UNB filter transmission curve and f\ is the 
flux density of the emission line based on the Rhoads et al. (2003) spectra. This takes the 
likely underestimation of Lya line flux into account. Following this, we were able to convert 
the convolved line flux into a magnitude by the use of the formula: 



and c is the speed of light. In order to ensure that all instrumental effects were taken 
into account, the last step in this process was to incorporate the detection fraction at each 
magnitude bin (Section I2.4( Figure |2]). This detection fraction was multiplied by the number 
of galaxies in each magnitude bin before we converted those magnitudes to Lya luminosities. 
The number of detected galaxies in each Lya luminosity bin was then used to estimate how 
many Lya emitters we should expect from our survey. 




(5) 



where 




(6) 



4 http : //www.noao . edu/ets/newf irm/documents/l056 , /o20nm / 20data°/ 20pack .xls 
http : //www.noao . edu/ets/newf irm/documents/1063y o 20nm / 20data°/ 20pack .xls 
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This Monte Carlo simulation was run for each individual NEWFIRM UNB filter and 
then iterated 10 times for each filter. We then averaged the results of those ten iterations, 
and the expected number of sources per filter at each magnitude bin is shown in Figure HI 
Looking at the limiting magnitudes for each filter and integrating out to that point, we find 
that we expect ~1 source per filter to be a true z = 7.7 Lya emitter. This result should be 
viewed with some caution, however, as we have assumed that the Lya luminosity function 
does not undergo significant evolution between redshifts of 6.5 and 7.7, and that the emission 
line profiles of all Lya emitters at z = 7.7 are the same as that of a Lya emission line profile 
at z — 5.7. Moreover, this simulation makes use of the Kashikawa et al. (2006) luminosity 
function, which is based upon detections in only one field, the Subaru Deep Field. We 
expect that there will be field-to-field variations among Lya emitters, and thus the expected 
number of sources may be different for the COSMOS field. Various methods can be taken to 
estimate the cosmic variance between fields of these Lya emitters; for large survey volumes 
(~2 x 10 5 Mpc 3 ), Tilvi et al. (2009) estimated this variance to be >30%. We can also utilize 
the method of Trenti & Stiavelli (2008) to calculate the cosmic variance expected assuming 
2 intrinsic Lya sources as reported by our Monte Carlo simulations. Using the redshift 
interval of 0.1 for our two filters and a Press-Schechter bias, the Trenti & Stiavelli cosmic 
variance calculate,! returns an expected cosmic variance of ~58% (this number drops by 
~10% with Sheth-Tormen bias). In add ition to cosmic yariance, we note here that the most 



recen t z=5.7 and 6.5 Lya LF results ( IHu et al.l I201CH : lOuchi et al.l l2010t Kashikawa et al. 



20111 ) indicate that the predictions for the number densities of z=7.7 emitters based on the 



results of Kashikawa et al. (2006) may be up to a factor of 3 too large. The uncertainty 
that remains in the lower- z Lya LFs, taken in concert with the significant role that cosmic 
variance can play, requires the Monte Carlo simulations to be viewed with caution rather 
than interpreted at face value. 



6. z = 7.7 Lya LUMINOSITY FUNCTION 

The primary goal of this work is to constrain the observed luminosity function (LF) of 
Lya emitters at z = 7.7. As described in the Introduction, our survey has the advantage 
of being both wide and deep, leading to our volume and limiting flux being better than or 
comparable to most other surveys at this redshift or higher (Tabled]). To construct our z 
= 7.7 Lya LF, we have used our four targets, which have Lya line fluxes of 12.1, 8.6, 8.1, 
and 8.0 x 10~ 18 erg s _1 cm -2 , all higher than our survey's limiting flux of 8 x 10~ 18 erg 
s _1 cm -2 . Using our survey volume of ~1.4 x 10 4 Mpc 3 per filter, we have constructed a 



http : / / casa. Colorado . edu/^trenti/CosmicVariance .html 
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cumulative Lya LF which is shown in Figure |5j In order to calculate the errors, we have 
used Poissonian statistics with Bayesian likelihood (where the likelihood is equivalent to the 
probability), assuming each target is in a separate luminosity bin. The assumption in this 
error model is that the luminosity distribution is uncorrelated - one galaxy having a certain 
luminosity has no effect on the luminosity of another - and thus it is safe to divide our 
luminosity bins this finely. As we have a prior probability distribution that is uniform in 
the expected number m (from Poisson), then the lower and upper limits on our errors are 
both finite. We have plotted our LF along with the results of several other major surveys. 
Two of these, Tilvi et al. (2010) and Hibon et al. (2010), focus on the same redshift as 
our work but have made use of alternate fields (CETUS and CFHT-LS Dl, respectively) 
and thus cosmic variance may come into play. Neither survey has published spectroscopic 
confirmation at this time, although Clement et al. (2011) note that VLT spectroscopy failed 
to detect Lya emission from the five most luminous Hibon et al. (2010) candidates, a result 
which will be published soon; we plot the Clement et al. (2011) upper limits for z = 7.7 
Lya LFs as well. We also plot the spectroscopically confirmed lye et al. (2006) z = 6.96 
Lya emitter. In addition to individual data points, we also plot three curves from previous, 
slightly lower redshift surveys, which are based on best-fit functions to spectroscopically 
confirmed data. These include the Ouchi et al. (2008) z = 5.7 LF, as well as the well-cited 
Kashikawa et al. (2006) z = 6.5 LF data. We have also included a new data set, the z = 
6.6 LF from Ouchi et al. (2010), which samples a different field from the Kashikawa z = 6.5 
LF but agrees well with that data set. 

Our data points agree quite well with the findings of Tilvi et al. (2010), despite our 
surveys probing different fields, but only the most luminous target matches with the Clement 
et al. (2011) upper limits for z = 7.7 within the error; our data points also may seem 



incon sistent with recent Lyman break surveys (e.g., iPentericci et al. 



2011; Schenker et al. 



20111 ). Assuming that all four of our data points are real Lya emitters at z = 7.7, we find 
good agreement with the z = 5.7 Ouchi et a l. (2008) function, yet our fainter luminosity 
objects do not ag ree with the z = 6.5 LF faashikawa et al" 2006 ) nor the z = 6.6 LF 



(lOuchi et al.ll2010l ). This could imply moderate evolution between z = 7.7 and z = 6.6, yet 
it is interesting that in such a case, the LFs at z — 7.7 and z = 5.7 would be in agreement. 
If we instead take the more conservative approach indicated by our Monte Carlo simulation 
results and assume only the two most luminous of our targets to be real Lya emitters, then 
our results are consistent within error with both the z = 5.7 Ouchi et al. (2008) function and 
the z = 6.5 Kashikawa et al. (2006) function. Given that it is more probable statistically 
that only these two candidates are real, we must take this case as the more likely outcome 
until we have the ability to spectroscopically verify our sources. Thus, given the agreement of 
our two most luminous emitters with the LFs at both lower redshifts, there is no conclusive 
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evidence for evolution of the Lya LF over the redshift range 5.7 < z < 7.7. 

7. SUMMARY 

We have utilized two custom-made UNB filters, at wavelengths 1.056 and 1.063 /im 
(FWHM 7.4 and 8.1 A, respectively) on the NEWFIRM camera at the KPNO 4m Mayall 
telescope to perform a deep and wide search for z = 7.7 Lya emitters in the COSMOS field. 
Our study comprised a co-moving volume of 2.8 x 10 4 Mpc 3 (survey area ~760 arcmin 2 ) 
and probed down to a limiting flux of ~8 x 10 -18 erg s _1 cm -2 (50% completeness limit), 
which is comparable to or better than previous Lya searches at similar redshifts. 

We used a very detailed selection procedure, making use of five different quantitative 
parameters as well as qualitative methods to ensure narrow-line detection and elimination 
of contaminants (including 3 low- 2 interlopers already known in the COSMOS catalog). We 
were left with a total of four candidates (three detected in the 1.056 /im filter, one in the 
1.063 jj,m filter), each detected at 5 sigmas or higher, down to line fluxes of 8 x 10 -18 erg s _1 
cm -2 . Detailed Monte Carlo simulations would suggest that up to two of these candidates 
are real. If we assume only two real candidates, comparison of the resultant Lya LF to the z 
= 5.7 Ouchi et al. (2008) and z = 6.5 Kashikawa et al. (2006) Lya LFs would indicate that 
there has been no significant evolution of the Lya LF between 5.7 < z < 7.7. This result 
is consistent with the findings of Tilvi et al. (2010) - that work used the same instrument, 
equivalent field of view, reached similar flux limits, but was half the volume of our work, 
owing to our extra filter. 

To pin down the neutral hydrogen fraction at z = 7.7 and thus the stage of the reion- 
ization process at that epoch, we will need more detailed Lya LFs in order to accurately de- 
termine the characteristic luminosity, L*, of these objects. This will require a) spectroscopic 
confirmations of these candidate high-redshift Lya emitters, and b) surveys encompassing 
more fields, as cosmic variance is likely to affect the number of emitters found in each region. 
The success of previous narrowband surveys at identifying Lya emitters at lower redshift and 
the robustness of our current set of candidates mean that the future of such high- z studies 
is quite promising. It is also encouraging to note that our results match well with surveys 
in different fields. The pursuit of surveys such as these, along with the advances that should 
be brought to the field by JWST and other new instruments, should provide a bright future 
for the study of the Dark Ages. 
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Table 1. Comparison to previous high- z Lya searches 



z 


Survey Vol. (Mpc 3 ) 


Detec. Limits (erg s 1 cm 2 ) 


# of LAE Detec. 


Refs 


(1) 


(2) 


(3) 


(4) 


(5) 


7.7 


6.3xl0 4 


8.3xl0~ 18 


7 


Hibon+ 2010 


7.7 


1.4xl0 4 


7xl0" 18 


4 


Tilvi+ 2010 


7.7 


2.8xl0 4 


8X10" 18 


4 


This work 


8-10 


35 


2X10" 17 


6 


Stark+2007 


8.8 


6.3xl0 4 


1.3X10" 17 





Cuby+2007 


8.96 


1.12xl0 6 


6xlCT 17 





Sobral+2009 



Note. — Col.(l): Lya redshift probed. Col. (2): Survey volume. The volume of our study is 
2x deeper than Tilvi et al. 2010 owing to our use of a second filter for candidate selection (each 
filter probes a volume of 1.4xl0 4 Mpc 3 ). Col. (3): Survey flux detection limit. Col. (4): Number of 
candidate Lya emitters detected. Col. (5): Survey reference. 



Table 2. Properties of four candidate Lya emitters. 



Cand. # RA Dec MagAB Line Flux (erg s 1 cm 2 ) L Lya (erg s- 1 ) EW LyQ (A) 



(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


1 


10:00:46.94 


+02:08:48.84 


21.87 


1.21x10- 


17 


8.34xl0 42 


7.32 


2 


10:00:20.52 


+02:18:50.04 


22.25 


8.55x10" 


18 


5.90xl0 42 


5.17 


3 (a) 


09:59:56.21 


+02:10:09.84 


22.31 


8.09x10- 


18 


5.58xl0 42 


4.89 


4 


10:00:48.79 


+02:09:21.24 


22.39 


8.00x10- 


18 


5.52 xlO 42 


4.84 



Note. — (a) Candidate #3 fails the optical selection criterion following an aperture correction of the Subaru data; 
it is thus considered more uncertain than the others. Col.(2): RA in J2000. Col. (3): Dec in J2000. Col. (4): AB 
magnitude, calculated using isophotal flux in SExtractor. Col. (5): Line flux, in the 1.056 /iia band for candidates 
1-3 and in the 1.063 /im band for candidate 4. Col. (6): Lya luminosity, calculated assuming a redshift of z=7.7. 
Col. (7): Lower limit rest- frame equivalent width of Lya emission, calculated as described in Section 13.31 equivalent 
widths. 
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Lya Redshift 

7.60 7.65 7.70 7.75 7.80 7.85 
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Wavelength (fim) 



1 00 1 ^ ^ I — 



Fig. 1. — NEWFIRM ultra-narrowband filter profiles, centered at 1.056 and 1.063 pm. The 
solid red, dashed, and dotted lines show the range of profile from center (solid red) to ~ 40% 
(dashed) and 80% center-to-corner distance (dotted) of the field of view. The Rousselot et 
al. (2000) OH sky spectrum is also plotted. 
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21.0 22.0 23.0 24.0 21.0 22.0 23.0 24.0 

1.056 Magnitude 1.063 Magnitude 

Fig. 2. — SExtractor detection percentages for randomly generated false sources in 0.1 mag- 
nitude bins in a given UNB band. The light gray bar marks 50% completeness. That value 
is used as the limiting magnitude (AB) for our survey. Our 1.063 /im band data is 0.1 mag 
deeper than our 1.056 ^m band data. 
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Fig. 3. — Image cutouts for our four Lya emitter candidates. UNB cutout images have 
been Gaussian smoothed according to the average seeing of the 1.056 and 1.063 [xm data. 
Cutouts are 50" on each side, and the circles are 16" in diameter (corresponding to ~400 kpc 
at z = 7.7). Each row represents one candidate. The optical column shows a weighted chi- 
squared combination of B, r, and i band images from Subaru. The middle two columns show 
our UNB NEWFIRM data. The J band column represents data from UKIRT. Candidates 
#1-3 are detected in 1.056 /im but not in any other band; candidate #4 is only detected in 
1.063 /im. All other objects present only in one UNB band are either transients (detected 
in only one yearly UNB stack) or fail to meet the Lya selection criteria. 
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Fig. 4. — Average results of the ten Monte Carlo simulation iterations, showing the number 
of Lya emitter sources we should expect at each magnitude bin. Integrating out to our 
limiting magnitudes for each field (Section 12.4ft . indicated by the vertical line, we expect 
roughly 1 source per filter to be a real z = 7.7 Lya emitter if there is no evolution in the 
luminosity function between z = 6.5 and z = 7.7. 
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Fig. 5. — Lya luminosity function compiled from data from this work (dark blue squares) and 
those previously published. Errors are calculated us ing Poissonian statis tics with Bayesian 
likelihood. The light blue triangles (jTilvi et al.ll2010[ ) and orange circles ( IHibon et al.ll2010f ) 
represent candidates at z=7.7 that have not been spectroscopically confirmed. Recent follow- 
up spectroscopy by Clement et al. (2011) did not detect Lya emission from the five most 
luminous candidates of Hibon et al. (2010), so we have excluded these objects from the 
present figure. The black dotted curve represe nts upper limit s on the z=7.7 luminosity func- 
tion from Clement et al. (2011). The red X ( live et al.l 120061 ) represents a spectroscopically 
confirmed z=6.96 Lya emitter. The three black dashed and dashed-dotted curves represent 
best fit luminosity functions to spectroscopically confirmed Lya emitters at redshifts z=5.7 - 
6.6. Note that our Lya luminosity function agrees well with that of Tilvi et al. (2011). Our 
two most luminous candidates are consistent within the error with either the z=5.7 or z=6.5 
luminosity function, and thus we do not see any evidence for significant evolution of the 
Lya luminosity function or neutral hydrogen fraction of the IGM between z=5.7 and z—7.7. 



